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ATMOSPHERIC  DISPERSION  RELATED 
TO  THERMAL  SPOT  TRACKING 

Lu  Qisheng,  Liu  Zejin,  Jiang  Zhiping, 

Zhang  Zhengwen,  and  Zhao  Yijun 

Department  of  Applied  Physics 

National  University  of  Defense  Technology 

Changsha,  Hunan  410073 

ABSTRACT  The  propagating  traces  of  two  light  beams  of  A,  =  3.8/rra  and  I0.6/jm  were  calcu¬ 
lated  in  the  atmosphere  with  spherical  symmetrical  index  distribution  of  refraction,  and  the  influ¬ 
ence  of  the  chromatic  dispersion  on  the  thermal  spot  tracking  was  discussed  in  this  paper. 

KEY  WORDS  atmospheric  dispersion,  light  trace,  light  spot  lock-in,  thermal  spot  tracking. 

I.  Introduction 

As  light  propagates  in  a  medium,  the  refractivity 
distribution  of  the  medium  strongly  affects  the  propagation  path 
of  light.  Refractivity  is  related  to  the  factors  of  the  density, 
temperature  and  pressure  of  the  medium.  Moreover,  refractivity 
is  also  related  to  factors  of  terrestrial  gravitation,  space 
environment,  meteorology,  terrain,  and  ecological  environment. 
There  have  been  numerous  research  achievements  about  these 
effects,  compiled  in  related  publications  [1,2].  If  only  the 
terrestrial  gravitation  is  considered,  the  distribution  of 
refractivity  with  atmospheric  height  can  be  assumed  as 
spherically  symmetrical.  Especially  in  some  local  environments. 
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this  approximation  is  basically  established. 

At  present,  the  problems  of  aim.ing  and  tracking  with  laser 
beam.s  are  often  involved  in  scientific  research  and  defense 
technology.  In  particular,  the  precise  tracking  of  m.oving 
objects,  and  steadily  locking  with  laser  light  spots  onto  a 
certain  site  of  a  m.oving  object  is  very  im.portant.  This  precise 
tracking  should  adopt  a  closed  cycle  system,  to  tim.ely  revise  the 
tracking  error  in  order  to  satisfy  high-precision  tracking 
requirements.  When  a  moving  object  is  illum.inated  with  a  laser 
beam,  two  kinds  of  radiation  will  be  generated:  one  is  diffused 
reflected  light,  and  the  other  is  the  therm.al  radiation  due  to 
heating  by  the  laser.  Both  kinds  of  radiation  can  be  considered 
as  the  feedback  signal  of  a  closed-ring  chain.  The  feedback 
signal  using  diffused  reflected  light  is  called  echo  tracking. 
The  feedback  signal  using  therm.al  radiation  is  called  therm.al- 
spot  tracking.  The  basis  for  selecting  the  tracking  schem.e  is 
how  to  upgrade  the  signal-to-noise  ratio  and  to  shorten  the 
response  tim.e,  am.ong  other  elem.ents.  For  example,  the  response 
is  slow  in  the  case  of  therm.al-spot  tracking  as  tens  of  m.s  are 
required  to  build  up  therm.al  spots.  In  the  case  of  echo  waves, 
the  transm.ission  tim.e  of  light  in  air  is  the  only  consideration. 
In  a  shared-aperture  tracking  system,  the  diffused  reflection  of 
the  system,  becomes  the  sam.e  frequency  interference  for  echo 
tracking.  However,  in  the  case  of  thermal-spot  tracking,  the 
wavelength  of  therm.al  radiation  can  be  used  as  the  light 
V7avelength  of  the  feedback  signal.  Just  by  adopting  the  m.easure 
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of  light  filtering,  a  relatively  high  signal-to-noise  ratio  can 
be  obtained. 

In  the  scheme  of  thermal-spot  tracking,  the  wavelength  of 
the  emitted  light  differs  from  the  wavelength  of  the  received 
signal  light.  With  the  effect  of  chromatic  dispersion  in  the 
atmosphere,  these  two  kinds  of  light  propagate  along  the 
different  paths;  this  will  affect  precise  tracking.  In  the 
following,  l.Spm  is  used  as  the  wavelength  of  the  emitted  light, 
and  10.6pm  is  used  as  the  wavelength  of  the  received  light. 

Thus,  an  example  is  used  to  calculate  the  deviation  of  the 
transmission  locus  in  the  atmosphere.  The  effect  of  such 
deviation  on  thermal-spot  tracking  is  discussed. 
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II.  Propagation  of  Light  Rays  in  Spherically  Symmetric  Atmosphere 
Assume  the  refractivity  distribution  in  the  earth's 
atmosphere  is  related  to  height  h  above  the  ground,  that  is 

n  (/•)  =n  (r,+  A)  (1) 

In  the  equation,  r^  is  the  earth's  radius.  When  a  laser  beam  is 
emitted  upvrard  from  A  on  the  ground  with  angle  <Dj  (refer  to 
Fig.  1)  ,  its  path  can  be  calculated  vyith  the  light  ray  equation 


[3]  : 


f'  dr 


(2) 


In  the  equation 
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Fig.  I  The  geometry  of  light  line 

propagation  in  the  atmosphere 


(3) 

«,(r)  =  l  +  A^,(r)  (4) 

AT  (0.7748  + 0.0046976A-')x  10 (5) 

"  t(r) 

Eq.  (5)  is  adapted  to  the  environment  [1,2] 

;7(r,)=  1.01325  X  lO’Pa  ,  T (r,)  =  288.15K  .  In  the  equation,  the 

atmospheric  pressure  distribution  p(r)  with  height  has  Pa  as  the 
unit;  and  the  wavelength  is  given  in  pm.  To  calculate 
quantitatively,  in  Eq.  (5)  p(r)/t(r)  is  approxim.ate  with  the 
ideal  gas  state  equation,  that  is 

P  (r) 

In  the  equation,  k  is  Boltzm.ann's  constant. 

By  using  the  results  of  Appendix  2  in  [4] ,  an  analytical 
function  is  (by  trial  fitting) 
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P  =  P  (O  X  10-*'” 


(7) 


In  the  equation,  h  is  the  height  from  the  ground,  given  in  km>. 
The  condition  is  adaptable  only  within  the  range  h=0  to  120km.. 
Since 

p(0/T  (0  =  351.64 


therefore 


finally,  we  obtain 


P(r+h) 

T{r,+  h) 


«351.64x  10-*'” 


(8) 


«i(^+/i)  =  l  +  (272.45+1.652/l-')  x  10-<‘*4> 


(9) 


Substitute  Eqs .  (3)  and  (9)  in  Eq.  (9)  to  obtain  the  light 

path.  It  can  be  proved  that  the  light  path  is  a  planar  curve  [3] 
in  the  plane  included  with  angle 


III.  Calculation  Results  and  Discussion 

Assum.e  Aj=3.8tJin  and  A,2=10 . 6|jm. ,  and  the  initial  em.ission 
angle  is  45",  rg=6000km..  From,  calculations,  the  laser  (of  two 
wavelengths)  locus  in  air  is  shown  in  Fig.  2. 

In  the  figure,  the  abscissa  A",. =  (/•,+  /!) sin(?>.,.  ,  and  the  ordinate 
^li- h)cos6^^-r^.  and  indicate  the  projections  on  the  x- 

and  y-axes.  A=3.8|jni  and  10.6pm..  The  laser  emits  a  locus  with 
<J>g-45".  Fig.  2  exhibits  and  .  These  are 

expressed  in  terms  of  deviations  from,  the  actual  locus  when 
em.itted  at  45'^.  This  is  to  am.plify  the  details  of  the  deviation. 
Although  in  Fig.  2,  m  is  used  as  the  unit  of  the  vertical  axis, 
we  are  still  unable  to  separate  the  light  locus  of  the  two 
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Fig.  2  The  propagating  trace  of  / ,=  3.8Mm  and  Fig.  3  The  difference  between  of  the  propagating  trace 
^  of  the  A,  =  3.8#im  and  10.6/tm 


wavelengths.  Therefore,  we  plotted  Fig.  3.  In  Fig.  3,  the 
vertical  coordinate  has  inin  as  the  unit  in  terms  of  _  y.. 

In  addition,  the  numerical  results  corresponding  to  h=l,  3,  10, 
100,  and  120km,  are  listed  in  Table  1. 

The  propagation  path  of  light  is  reversible  in  a  medium..  If 
from,  point  A,  a  laser  emits  at  the  vfavelength  3.8(jm.  toward  point 
B  at  the  angle  O,,  to  burn  a  therm.al  spot  at  point  A,  the  10.6|jm. 
radiation  emitted  by  the  therm.al  spot  returns  to  point  B  with 
different  light  paths.  As  revealed  from,  the  calculation  results 
m.entioned  above,  the  incident  angle  toward  point  A  in  the  path 
of  10.6|jm  irradiated  light  is  smaller  than  If  the  system,  is 

so  controlled  that  the  laser  {3.8|Jin)  emits  at  the  angle  <t>\,  then 
the  laser  light  spot  will  m.ove  upward  from,  point  B.  In  such 
repeated  cycles,  the  light  spots  cannot  be  locked,  which  will 
autom.atically  leave  the  target.  Angular  revision  should  be  m.ade 
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Table  1  The  data  of  light  line  trace  with  A|  =  3.8Mm  ,  lO.hfim  .  d)  =  45“ 


hj  km 

A'„/km 

mm 

1 

0.99995 

-3.532  X  I0“- 

0.99995 

-3.534  X  10-= 

-  2.372  X  10“= 

3 

2.99954 

-2.913  X  10“' 

2.99954 

-2.914X  IQ-' 

-3.364  X  10’’ 

10 

9.99415 

-2.461 

9.99415 

-2.461 

5.338  X  10“' 

100 

99.24261 

-49.962 

99.24259 

-49.948 

14.636 

120 

1  18.90544 

-60.654 

1  18.90542 

-60.636 

17.826 

if  precise  locking  is  intended.  If  the  target  object  moves  at 
the  velocity  v,  then  the  tracking  system  should  rotate  the  light 
beam  at  the  angular  velocity  iv/R!  (with  the  direction  Rxv) . 
Moreover,  revisions  are  made  from  tim.e  to  tim>e  on  the  path  (or 
the  emission  angle) .  Only  Eg.  2  is  subjected  to  integration  over 
tim,e,  and  we  can  obtain  d©/dt,  as  well  as  the  relationship 
between  the  velocity  v  and  d<Dj/dt.  When  adopting  the  precise 
tracking  schem.e  of  therm.al  spot,  we  should  establish  the 
relationship  between  dO'ydt  and  dO^/dt.  This  is  not  the  problem, 
to  be  solved  in  this  article,  therefore  this  is  omitted. 

The  first  draft  of  the  article  was  received  on  April  2, 

1992;  the  final  revised  draft  was  received  for  publication  on 
June  11,  1992. 
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FOLDING  FFT  MODIFICATION  FOR  LASER  PROPAGATION 
IN  THE  TURBULENT  ATMOSPHERE 

Wang  Yingjian  and  Wu  Yi 

Anhui  Institute  of  Optics  and  Fine  Mechanics, 

P.O.  Box  1125 
Hefei  230031 

ABSTRACT  The  traditional  method  of  numerical  simulations  for  laser  beam  propagation  in  the 
turbulent  atmosphere  is  modified.  A  folding  FFT  treatment  is  put  forward.  High  spatial  frequency 
terms  of  the  turbulent  spectrum  are  sampled  and  low  spatial  frequency  terms  are  revised  using 
this  treatment.  Numerical  simulation  results  show  that  the  calculation  speed  and  accuracy  are  both 
enhanced  and  the  required  computer  memory  is  reduced.  Also  the  numerical  results  reflect  the  real 
properties  of  strong  and  weak  effects  of  the  atmospheric  turbulence. 

KEY  WORDS  atmospheric  turbulence,  laser  propagation.  FFT  method. 

I.  Introduction 

Numerical  calculations  of  the  effects  of  turbulence  in  the 
atmosphere  are  one  of  the  important  aspects  in  studying  the 
atm.ospheric  transmission  of  lasers.  These  calculations  relate  to 
establishing  a  four-dimensional  calculation  program  of  the  high- 
powered  laser  transm.ission  in  the  atmosphere.  The  conventional 
sim.ulation  algorithm,  of  the  turbulent  phase  screen  in  the 
atmosphere  is  the  lim.itation  by  the  computer  computational  speed 
and  internal  memory.  It  is  very  difficult  to  ensure  the 
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requirements  of  computational  accuracy  with  regard  to  the 
turbulence  spectrum  as  well  as  the  total  width  of  space  sampling 
and  sample  density.  Therefore,  the  information  lost  in  the  high- 
and  low-frequency  spectra  of  turbulence  is  very  great.  With 
regard  to  the  optical  effects  of  such  loss,  this  involves  the 
tilt  of  the  light  beam  and  scintillation  of  light  intensity.  It 
is  even  so  with  the  intensive  turbulence  effects.  Therefore  some 
authors  calculated  turbulence  effects  by  using  various  modified 
spectra  [1,2].  The  authors  improved  on  the  conventional 
simulation  algorithm  of  the  turbulence  phase  screen  in  the 
atmosphere  and  proposed  a  folding  FFT  algorithm,  thereby  solving 
the  problem  of  numerically  describing  the  turbulence  spectrum  in 
the  calculation  method.  As  revealed  in  the  results  of  the 
numerical  simulation  of  laser  beam  transmission  in  a  turbulent 
atmosphere,  this  algorithm  upgrades  the  computational  speed  and 
accuracy,  lowers  the  requirements  on  computer  internal  memory, 
and  more  truly  simulates  the  strong  and  weak  effects  of 
atmospheric  turbulence.  In  the  authors'  view,  this  approach  is 
helpful  to  constructing  a  four-dimensional  computer  program  of 
high-powered  laser  transmission  in  the  atmosphere,  especially 
when  considering  the  interaction  of  turbulent  flow  and  thermal 
blooming . 

II.  Inversion  of  Turbulence  Phase  Screen 

Generally,  numerical  simulation  on  the  effects  of 
atmospheric  turbulence  in  laser  transmission  is  carried  out  by 
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using  the  phase  screen  approximation  [1-4].  In  this  approach, 
the  turbulent  atmosphere  is  divided  into  several  layers  along  the 
direction  in  winch  the  light  beam  is  transmitted.  The 
contribution  of  each  layer  amounts  to  ciianging  the  phase  of  the 
transmission  light  beam  in  that  layer;  however,  the  light  beam  is 
considered  as  transmitted  in  vacuo  when  the  beam  is  between  two 
layers.  For  the  atmospheric  turbulence  layer  with  thickness  Ciz , 
the  action  on  the  phase  perturbation  difference  of  the  light  beam 
is: 

Acp^k  5n(r,  z')dz'  /i  n 


In  the  equation,  k=2n/A  is  the  wave  number;  5n  is  the 
refractivity  perturbation  of  the  turbulent  atmosphere;  z  is  the 
transmission  direction;  r  is  the  coordinate  of  the  point  as  the 
light  beam  passes  through  the  cross-sectional  plane.  The 
approximate  conditions  of  the  phase  screen  are  as  follows  [1-4] : 
Dz  is  longer  than  turbulence  correlation  distance,  but  smaller 
than  rgVA;  r^  is  the  transverse-direction  coherence  length  of  the 
turbulence.  Due  to  the  randomness  of  atmospheric  turbulence,  the 
phase  screen  is  unable  to  be  executed  by  using  Eq.  (1).  The 
authors  adopt  the  phase  correlation  function  and  density  of  the 
refractivity  perturbation  spectrum  to  invert  the  phase  screen. 

The  self-correlation  definition  of  the  phase  screen  is  [5] : 
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>■,)  =  <  A(P  (r,  ,  r)A(p(r,,  z)  > 


-  <  (^«(r, ,  r,  )(5«  (*■;,  r, )  >  Jr, 

J  *}  - 


<  >  indicates  the  overall  average  of  the  system.  Based  on  [1] 
and  [5] ,  there  is ; 

<5n  (r,,  z,  )dn{r„  z, )  >  =  It:  C;  (z, ,  z,)d  iz,  ,  z,)  (^(k,  =  0,  T) 

J  J  ~X 

■  exp  (iK  ■p)d-T  (3  ) 

C^j,  (z)  is  the  vertical  distribution  of  the  structure  constant  of 
turbulence  ref ractivity ;  5(z)  is  the  Dirac  function; 

^=r| -r„  iJc^l  =  (k-  + -  is  the  conjugate  variable  of  the  Fourier 
transform  domain  of  p;  o,(k.  =  o, /T)  is 

the  density  function  of  the  refractivity  fluctuation  spectrum. 
With  respect  to  the  Kolmogorov  spectrum  [6] : 


O,  (k)  =0.033x-  5 


2n  5.92 


Ig  and  Lg  are,  respectively,  the  internal  and  external  dimensions 
of  the  turbulence.  With  respect  to  Von  Karman  spectrum  [6]: 

( K  )=  0.033  47:7' Uh  i 


Substitute  Eg.  (3)  into  Eg.  (2)  and  we  obtain; 


-2xk-  =  0,  x- )exp  (/ Cl(z')dz 
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Definition : 


F.^I'k)  =27tA--0  Jk:.  ,  iT)  Cl(z')dz' 


(5b) 


We  can  see  that  and  are  the  conjugate  Fourier 

transform  functions.  Thus,  the  light  beam  phase  fluctuation 
screen  induced  by  the  turbulent  atmospheric  layer  with  a 
thickness  (5z)  can  be  generated  by  the  conventional  method  [1,2] 


Aw  (r,  z)  =  q 


g(K  )F'^  -  (K)e\p  (iK  ■  r)d- k 


(6) 


In  the  equation,  g(K)  is  the  complex  type  two-dimensional 
Gaussian  random  white  noise;  and  q  is  the  calibration  vector. 
With  respect  to  isotropic  turbulence,  and  aO  (r)  rely 

only  on  and  jrl  .  Therefore,  only  if  the  density  function 

Ojj  of  the  refractivity  fluctuation  spectrum  is  given,  the  phase 
fluctuation  ti<D(r)  can  be  obtained  by  using  Eq.  (6).  However,  we 
know  that  the  turbulence  spectrum  is  very  wide  in  scope,  within 
the  range  2n/h^  and  2n/l(j,  the  turbulence  seriously  affects  laser 
transmission.  If  the  Cauchy  2  fast  Fourier  transform  algorithm 
is  used  for  numerical  calculation  on  Eq.  (6),  in  order  to  attain 
the  required  calculated  sampling  density  (CiK=2n/h^)  or  ^=ljj ; 
usually,  Ig  is  in  millimeters,  and  Lq  is  approximately  10  or  more 
meters,  the  divergent  dot-matrix  number  of  the  phase  screen  will 
attain  8192x8192.  However,  the  numerical  calculated  dot  matrix 
number  in  the  transmission  equation  is  only  a  small  part  of  them 
We  can  see  that  this  is  very  uneconomical  for  the  traditional 
divergence  calculation  method  corresponding  to  the  space  and 
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Fourier  spectrum  domain  dot-matrix  number.  Moreover,  it  is  very 
difficult  to  have  this  much  computer  memory  and  computational 
volume.  Thus,  many  authors  sacrifice  the  accuracy  (reducing  the 
dot  matrix  number)  to  upgrade  the  computing  speed  and  reduce  the 
internal  memory  [2,3]. 

III.  Folding  FFT  Algorithm 

We  cite  an  example  with  twice-folding  to  discuss  the 
algorithm  in  Eq.  (6).  If  the  number  of  space  divergent  points  is 
2Nx2N  and  the  sample  interval  is  dx  then  the  divergent  points  of 
the  Fourier  spectral  domain  are  2Nx2N,  ^=2n/N^.  After 
divergence  of  Eq.  (9),  we  have: 

N-\ 

A(pim,n)='£  X  f  im',n')exp[i2itinn'  +  mm')/  N]  (7) 

n  '  *— ,V  m'^  -N 

In  the  equation,  f  ^  g{m' ,n' ,n' ),  L  =  NAx  is  the 

overall  width  of  the  space  sample.  The  scalar  factor  q=Chc’^. 

Let  us  rewrite  Eq.  (7)  to  become: 

N-\ 

A(p{m,n)=Yj  exp  ii2nnn' I  N)[Yj  f{m(,  n' )exp  {i  2n  m{ml  N ) 

n'  <=  —  N  m'w  — A? 

^-1 

+  E  (8a) 

m'  =  0 


Let  m2'=m'-N,  and  applying  the  periodicity  of  the  trigonometric 
function,  we  obtain 

.V-I  If -I 

A(p(m,  «)=E  expiiln  mm' !  N){  E  [f  {m'-N ,  n' ) +f  {m' ,n')] 

m  '  =  0.  m' 

■  Qxp  {i  2nm' j  N)}  (8b) 

Similarly,  by  conducting  a  transform  on  the  summation  equation  of 
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n',  finally  we  obtain: 


A  (p  {m,  n)  =  Yj  X  f  {m'  —  N ,  n'  N)  +f  {m'—  N,  n  ) 

m  '  =  0  n  '=  0 

+f{m\n'-N)+f  (/n'./i')]exp  [i2n(m'  m  +  n'  nj  N]  (g^) 

Obviously,  only  NxN  fast  Fourier  transforms  are  required  for 
computing  Eq.  {8c);  the  maximum  sampling  frequency  is  2r\/Cix.  We 
compare  the  folding  FFT  algorithm  with  the  conventional  algorithm 
and  discover  the  following:  by  one-to-one  correspondence  of  the 
number  of  points  in  the  space  and  frequency  domains  in  the 
conventional  algorithm,  if  the  space  frequency  is  of  the  same 
height  for  sampling,  the  dot  matrix  number  should  be  N2xN2=2Nx2N. 
In  other  words,  2Nx2N  fast  Fourier  transforms  should  be 
conducted,  and  the  number  of  multiplications  is  [7] : 

Y  (r+ 1 )  =2Af^  (/■+ 1 ),  A^=2'  y  is  a  positive  integer.  The  number 
of  summations  is  (r+1 )  =4N^ (r+1 )  .  The  number  of  multiplications 
in  the  folding  algorithm  is  {1/2)nV  and  the  number  of  summations 
IS  N  (r+4) .  Thus  we  know  the  number  of  multiplications  in  the 
folding  FFT  algorithm  is  one-quarter  that  of  the  conventional 
algorithm,  and  the  number  of  summations  is  also  greatly  reduced. 
Correspondingly,  the  computer  internal  memory  required  is  only 
one-quarter  as  under  the  conventional  algorithm.  In  the 
conventional  algorithm,  if  the  NxN  dot  matrix  number  is  used,  the 
sampling  highest  space  frequency  is  2p/2Dx,  which  is  one-half  as 
under  the  folding  algorithm;  thus,  the  computational  accuracy 
will  be  greatly  reduced.  Besides,  we  can  see  in  the  above- 
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mentioned  operations,  in  principle,  more  than  twice-folding 
computations  can  be  conducted  to  further  upgrade  the  sampling 
accuracy  of  the  turbulence  density  spectrum.  In  short,  with  the 
same  sampling  frequency  (that  is,  the  same  computational 
accuracy)  the  folding  FFT  algorithm  is  better  than  the 
conventional  algorithm  with  respect  to  computing  speed  and  the 
internal  computer  memory  required. 

IV.  Numerical  Computational  Results  and  Discussion 

We  applied  the  twice-folding  FFT  algorithm  to  conduct 
numerical  simulation  on  the  equivalent-level  homogeneous 
atmospheric  transmission  of  a  collimated  gaussian  laser  beam;  the 
refractivity  fluctuation  spectrum  applies  the  Von  Karman 
function.  The  sampling  total  width  L=2.56m,  N=128,  L[|=62.8m, 
A=1.0ijm,  ^z=100m,  at  the  1/e  power  point  of  the  gaussian  light 
beam,  the  beam  waist  radius  ag=0.5m,  and  the  structure  constant 
C^jj  of  the  turbulence  refractivity  is  10'^^  to  10'^^ 

Corresponding  to  weak  and  strong  turbulence  features,  the  maximum 
transmission  distance  is  10km.  Fig.  la  is  an  example  of  phase 
screen.  We  can  see  that  there  are  obvious  indications  in  the 
phase  screen  from  small-scale  phase  fluctuation  to  large-scale 
phase  tilt.  Fig.  lb  shows  a  comparison  with  the  theoretical 
results  and  the  phase  screen  reconstruction  phase  structure 
function  D{r).  The  "+"  indicates  the  average  results  of  ten 
phase  screens;  the  solid  curve  shows  the  theoretical  results. 
Their  agreement  indicates  a  great  improvement  over  the  results  in 
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Fig  .  I  (a )  An  example  of  phase  screen  [  Horizontal  axis  drawn  with  2  points  interval  on  (.v.>')grid] 
(b)  Phase  structure  function,  r  distance  of  2  points  on  grid 


[1] .  This  indicates  that  the  structure  functions  in  the 
conventional  algorithm  gradually  decreases  and  approaches  zero 
with  the  greater  intervals  r.  The  results  are  quite  different 
from  the  theoretical  results  of 

Fig.  2  shows  an  example  of  intensity  distribution  on  a  light 
beam  cross  section  at  10km.  a  is  the  three-dimensional  display 
and  b  is  the  intensity  isopleth.  From  the  figure  we  can  see  very 
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Fig.  3a  is  the  normalized  variance  of  the  light  intensity 
scintillation  varying  with  transmission  distance.  The  "+" 
indicates  the  numerical  simulation  result;  the  solid  curve  shows 
the  theoretical  results.  In  the  weak  scintillation  region, 


L8 


Fig  .3  (a)  Normalized  intensity  deviation  square  vj  propagation  distance 

( b )  Probability  of  intensity  scintilation.  aj-  2.4 


is  proportional  to  )  [8]  ,  which  is 

consistent  with  the  theoretical  results.  When  the  light 
transmits  to  a  certain  distance,  rll/.z=\  is  the  turbulence 

perturbation  Rayleigh  distance,  and  light  intensity  simulation  is 

,  2 

at  saturation.  Then  decreases  with  a  certain  slope  with 

•  .  1 

increase  in  fij.  This  result  also  agrees  with  the  theoretical 
estimate  [8] .  Fig.  3b  shows  the  probability  distribution  of 
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< I >  is 


light  intensity  scintillation-  I  is  light  intensity; 
average  light  intensity;  and  the  statistical  region  radius  is  aj. 
The  "+"  indicates  the  numerical  simulation  result;  the  solid 
curve  is  the  distribution  of  positive  theoretical  logarithms  [9] . 
Both  are  also  quite  close.  In  the  situation  of  greater  I/<I>, 
the  numerical  sim.ulation  results  have  greater  divergence.  This 
is  because  the  greater  the  I/<I>,  the  sm>aller  is  the  probability. 
Thus,  the  statistical  num.ber  of  points  is  also  smaller.  With 
respect  to  the  numerical  simulation  of  light  intensity 
scintillation,  previously  no  research  reports  had  been  published. 
This  is  possibly  due  to  the  fact  of  low  sam.pling  density  and  a 
small  number  of  high-frequency  com.ponents  of  the  structure  phase 
screen,  thus  neglecting  the  effect  on  light  intensity 
scintillation  due  to  high-frequency  spectral  turbulence,  so  that 
it  is  difficult  to  obtain  better  results. 

Generally,  the  Strehl  ratio  is  adopted  in  evaluating  light 
beam,  quality  during  laser  beam,  transm.ission  in  a  turbulent 
atm.osphere.  Fig.  4  shows  the  variation  of  Strehl  ratio  with  D/r^ 
for  the  collim.ated  gaussian  light  beam,  transm.itted  in  a  turbulent 
atm.osphere.  We  use  D  =  2yJ~2  .  rg  varies  with  c’j,  and 

transm.ission  distance  z.  The  in  the  num.erical  simulation 

results  also  agrees  considerably  with  the  solid  curves  [10]. 

As  mentioned  above,  in  the  authors'  view,  this  com.putational 
m.ode  m.ore  truly,  overall,  reflects  the  effect  on  laser  beam 
transm.ission  by  atm.ospheric  turbulence  because  whatever  the  light 
intensity  scintillation  or  the  light-beam  quality  Strehl  ratio. 
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Fig  .  4  Variation  of  the  Strehl  ratio  with  Dj 


the  comprehensive  effect  is  imposed  on  turbulence  from  all 
frequency  ranges  (large  and  sm.all  scales)  . 

V.  Conclusions 

In  the  folding  FFT  algorithm  proposed  by  the  authors,  there 
is  an  improvem.ent  in  the  conventional  numerical  computation  of 
laser  transm.ission  in  a  turbulent  atm.osphere.  The  folding 
algorithm,  has  higher  com.putational  accuracy  and  speed,  and  lower 
requirem.ents  for  internal  com.puter  m.emory.  To  invert  a  phase 
screen  on  a  VAX-II,  7m.in  is  required  under  the  conventional 
algorithm.,  but  only  a  little  more  than  3min  with  the  folding  FFT 
algorithm..  This  will  have  a  great  effect  on  the  atmospheric 
transmission  of  a  high-pov7ered  laser  (especially  when  considering 
the  interaction  between  therm.al  bloom.ing  and  turbulence)  ,  and  the 
establishment  of  a  f our-dim.ensional  num.erical  com.putational 
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program.  Tx-yice-f olding  can  yield  satisfactory  results. 

Actually,  the  folding  FFT  algorithm  can  be  used  not  only  in 
numerical  sim.ulation  of  laser  transmission  in  the  atmosphere,  but 
also  it  is  useful  in  computing  the  number  of  sam>pling  points  in 
the  f requency-dom.ain  that  requires  a  m.uch  larger  num.ber  of 
sam.pling  points  in  the  space  dom.ain,  or  vice  versa. 

The  authors  express  their  gratitude  to  comrade  Wang  Ningquan 
because  of  his  help  in  preparing  the  figures  and  tables. 

The  first  draft  of  the  article  was  received  on  December  28, 
1991;  the  final  revised  draft  was  received  for  publication  on 
July  28,  1992. 
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ABSTRACT  The  purpose  of  this  monograph,  which  is  divided  into  two  parts,  is  to  give  a  gener¬ 
al  review  for  the  recent  advances  in  beam  transformation  optics,  including  some  results  obtained 
by  our  group.  Part  I  deals  with  the  diffraction  theory  in  the  time-space  domains  and  operator 
optical  methods,  as  well  as  their  applications  to  laser  optics  and  high  power  laser  technology. 
KEY  WORDS  beam  transformation,  high  power  laser  technology,  diffraction  integrals  in  the 
time  — space  domains,  operator  optical  methods. 


I.  Introduction 

In  the  research  scene  dealing  with  laser  beams  passing 
through  optical  elem.ents  and  system.s  (including  optical 
resonator)  with  the  transm.ission  transf orm.ation  rule  and  the 
control  of  light  beam,  quality,  light  beam,  transf  orm.ation  optics 
(also  referred  as  to  as  laser  optics)  has  seen  rapid  developm.ent 
of  late  owing  to  advances  in  laser  science  and  technology, 
especially  high-powered  laser  technology.  Some  new  research 
directions  have  appeared.  There  are  m.any  papers  and  research 
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reports  in  China  and  abroad  on  this  field.  In  reference  [1],  the 
authors  discussed  some  fundamental  analytical  m.ethods  and 
applications  on  beam,  transf orm.ation  optics.  In  this  paper,  by 
including  the  related  problem  of  high-powered  laser  technology, 
more  detailed  descriptions  are  given  on  the  advances  in  beam, 
transf  orm.ation  optics.  This  paper  is  divided  into  tv7o  parts. 

Part  I  discusses  diffraction  theory,  Ham.iltonian  optics,  and 
operator  m.ethods  in  the  tim.e— space  dom.ain.  Part  II  describes  the 
matrix  optical  m.ethod  and  Lie  algebraic  theory.  Finally,  a  brief 
analysis  is  made  on  the  correlated  problem.s  on  the  internal 
relationship  in  the  theory  and  high-powered  laser  technology. 

The  paper  em.phasizes  the  explanation  of  physical  concepts  and 
methods.  Some  application  examples  cited  also  serve  this 
purpose.  Further  understanding  about  details  and  actual 
applications  can  be  found  in  the  references. 


II.  Diffraction  Theory  in  the  Space-time  Dom.ain 

2.1.  Kirchhoff  and  Wickens-Fennell  diffraction  integral 

As  is  well-known,  the  fundam.ental  form.ula  describing  the 
scalar  diffraction  field  is  [2] 


1  +  cos  sc 


)dx, 


(2.1  ) 


In  the  equation,  Ui(x^,yi)  and  U2(x2,y2)  are,  respectively,  the 
field  oscillation  amplitudes  at  the  point  of  origin  and  the  field 
point;  p  is  the  spacing  between  the  point  of  origin  and  the  field 
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point.  (l+cosa)/2  is  called  the  tilt  factor;  a  is  the  angle  of 
intersection  between  the  normal  line  at  the  point  (xl,yl)  on  the 
S-plane  and  p;  X  and  k  are,  respectively,  the  wavelength  and  V7ave 
number.  The  integration  is  performed  on  the  S-plane.  Generally, 
Eq.  (2.1)  is  called  the  Kirchhoff  diffraction  formula. 

When  the  Fennell  approxim.ation  is  established,  in  other 
words,  the  distance  L  between  the  diffraction  plane  and  the 
observation  plane  is  much  greater  than  the  diffraction  aperture 
and  the  linearity  in  the  observation  region 

p^L.cosfl^I  (that  is,L»  I,  I  )(2.2) 

Eq.  (2.1)  is  converted  into  the  Wickens-Fennell  diffraction 


integral 


ik 


J/,(x,  ,y,)exp{  —  [  ]dx,dy\  (2.3) 


If  we  further  increase  L  so  that  the  following  equation  is 
satisfied 


(2.4) 


Then  Eq.  (2.3)  is  sim.plified  as  the  Fraunhofer  diffraction 
integral 


,;•.)=  (- 


IL 


-exp[;—  (J'c'  +  Ji)] 


T  ik 

..y,)exp[  -  —{x,x^+y^^ )']  dx,  dy,  (2.5) 


This  is  easier  than  Eq.  (2.3)  for  solving  by  computer.  With 
respect  to  the  diffraction  problem,  of  large  space  domain  in 
classical  optics,  Eq.  (2.3)  is  suitable;  the  equation  can  also  be 
used  in  studying  that  the  free  space  exists  betv7een  the 
diffraction  plane  and  the  observation  plane,  such  as  transmission 
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and  phase  control  problems  of  a  laser  beam  in  the  atm.osphere  when 
such  effects  as  the  following  are  neglected:  nonlinearity  effect, 
observation  and  linear  scattering.  For  a  long  time,  researchers 
thought  that  Eq.  (2.3)  was  universally  established;  however, 
advances  in  laser  technology  raise  a  challenge  to  this  viewpoint. 
The  related  research  m>ainly  includes  two  aspects.  First, 
recalling  the  original  Kirchhoff  formula,  study  the  gaussian 
laser  beam,  without  side  axes,  or  light  beam,  of  random,  wavefront 
passing  through  an  asymm.etrical  m.edium.,  in  addition  to 
diffraction  by  a  diaphragm,  of  arbitrary  shape  [3,4].  Second, 

Eq.  (2.3)  is  generalized  to  deal  with  such  problem.s,  including 
the  diffraction  problem,  of  the  space  dom.ain  of  a  com.plex  optical 
system.,  which  is  Collins'  form.ula  [5].  Thus,  further  expand  the 
equation  to  study  the  m.ore  general  tim.e-space  dom.ain  diffraction 
[6] .  In  the  following,  the  latter  problem  is  emphasized  in  the 
discussion . 


2.2.  Time- space  dom.ain  diffraction  integral 
Recently,  Kostenbauder  proved  [6]  that 
diffraction  integration  form.,  in  the  broad 
tim.e  variable)  (t,  f  =  frequency)  and  space 


U,i.x,,r,)=r] 


Ujix; )exp[ 


the  Wic)cens-Fennell 
sense  (including  the 
variable  (x,  6)  is 

<5  £  V.Xj  ]k,  i//,  (2.6) 


In  the  equation,  the 


beam,  transmission) 


m.atrix  (describing  the  light 
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the  matrix 


-V 

y 

t 

e 

(p 

V 


In  the  equation,  the  4x4  matrix  should  be 


expanded  into  a  6x6  matrix.  Eq.  (2.6)  should  also  be 
correspondingly  expanded. 


2.3.  Collins'  form.ula 

When  discussing  only  the  space  diffraction  problem,, 
Eq.  (2.6)  is  sim,plified  as 


UAx,)^ 


KB 


)exp[  - -^4-(,^x;-2x,.x,+  Z).vj)]  dx^  (2.10a) 
IB 


(5#0) 


vT 


i  kC  Y 

exp(- (.6  =  0) 


(2.10b) 


In  the  equations.  A,  B,  C,  D  are  elem.ents  in  the  following 
transformation  m.atrix 


(2.11  ) 


Eq.  (2.10)  is  the  well-l^nown  Collins'  formula,  which  can  be  used 
to  study  the  transform, ation  problem,  of  com,plex  optical  system,s  in 
which  a  laser  beam,  passes  through  a  transf orm.ation 

•  (  ®  I 

matrix  \  /  ,  and  also  can  be  used  to  the  motion  rule  [7]  of 

C  D 

laser  beam,  in  an  optical  resonator.  In  addition,  well-known 
interference  and  diffraction  phenom.ena  in  classical  optics  can 
also  be  unifiedly  processed  with  Eqs .  (2.10a)  and  (2.10b). 
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2.4.  Tims -domain  diffraction  integral 

Now  the  transmission  problem,  of  dispersion  optical  pulses  is 
studied.  Here  the  tim.e-space  domain  diffraction  integral, 

Eq.  (2.6),  is  transform.ed  into  the  tim.e  dom.ain  diffraction 
integral 


C/,(n)=  - 


InB 


t/i  (T,)exp[ (.4Tj-2T,T,+  Z)T;)]d'T,  ,(5#0)  (2.12a) 

J.£> 


V  A 


e'‘^'»exp( Ct;)  t/,(  ^  ),(5  =  0) 


(2.12b) 


In  the  equations 


T==/-^'(ajp)2=/-  -f  , 


(2.13) 


Vj  is  the  group  velocity;  6"  describes  the  dispersion  of  group 
velocity;  is  the  center  frequency  of  the  dispersion  pulse, 
however 

A  B  \  /  r, 


T, 

dx,  I  = 


C  D 


dx, 


(2.14) 


Obviously,  Eq.  (2.12)  can  be  readily  derived  from.  Eq.  (2.10) 
V7ith  the  following  tim.e-space  variable  analogy  [9,10]: 

^  .  t-p'z  =  x^ — ►  .r,  (2.15  ) 

More  rigorously,  proof  should  be  exercised,  beginning  from  the 
time  dom.ain  dispersion  wave  m.otion  equation  [11].  Take  note, 
similarly  as  in  the  classical  optics  situation,  for  sake  of 
convenience  in  studying  the  problem,  that  there  are  two  kinds  of 
equivalence  symbol  utilization  formulas  in  Eqs .  (2.7),  (2.8),  and 
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For  example,  another  equivalence  form  of  Eq.  (2.12a)  is 


.12) 


6',  (T,)exp[ iAx]-2x,r,+  Dtl)]dr,  ,(5#0)  (2.16) 


III.  Ham.iltonian  Optics 

Another  route  in  studying  light  beam  transformation  is  to 
try  to  include  the  m,ovem.ent  rule  of  the  light  frequency 
electromagnetic  field  into  the  framework  of  unified  theory  of 
physics.  Thus,  the  matured  physical  equations  and  the  related 
m.athem.atical  means  in  research  can  be  directly  used  to  deal  with 
problem.s  of  laser  optics.  Here,  analogy  and  inference  are 
m.ethods  in  scientific  research  with  another  successful  exam.ple. 
First,  let  us  recall  Ferm.at's  principle,  which  explains  that  the 
variation  is  zero  along  the  real  optical  path  of  the  light  rays. 
In  other  words. 


5 


f/j  {x,y,z)(i  1=0 


(3.1) 


In  the  equation,  n  is  the  refractivity  of  the  m.edium. .  Obviously, 
Eq.  (3.1)  is  sim.ilar  in  form,  to  the  Ham.iltonian  principle  in 
analytical  mechanics.  Next,  the  Helmholtz  equation  can  be 
directly  derived  from.  Maxwell’s  equations: 

U+k^  U=0  (2  2) 


with  side  axis 
become  [12] 


approxim.ation ,  the  equation  is 


dW  d^U 

dx-  dy- 


■2  ik 


dU 

2z 


transform.ed  to 

(3.3) 
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Eq.  (3.3)  is  generally  referred  to  as  the  equation  similar  to  the 
Schroedinger  wave  equation,  with  similar  forms  as  the 
Schroedinger  equation  well-known  in  quantum  m.echanics. 

Therefore,  physicists  well  acquainted  with  classical  m.echanics 
and  wave  motion  mechanics  will  very  easily  adopted  the  analogical 
approach  in  setting  up  Hamiltonian  optics  [13,  14]  and  operator 
optics,  with  elaboration.  In  Hamiltonian  optics,  there  are  the 
follov7ing  m.ost  im.portant  results: 

3.1.  Optical  regular  variable  and  operator  form 

The  positional  direction  of  a  light  beam,  on  a  transverse 
plane  can  be  written  as  the  following  equation  by  using  the 
regular  coordinates  q^^  and  q^  and  the  regular  m.om.enta  p^  and  p^ . 


In  the  side  axis  approxim.ation ,  there  is 

dx  dy 

=ne,p=n—  =ncp 

The  corresponding  regular  conjugate  operator  is 


(3.5) 


(3.6) 
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3.2.  Optical  Lagrangian  equation 


ei 


dz 


dz 


dL 


)  = 


)  = 


8L 

dx 


dL 


(3.7) 


8(  ^ 

dz 


In  the  equation,  the  optical  Lagrangian  argument  is 


L  =  n{x,y,z)  / 1  +  (  ^ 

V  az  rfz 


(3.8) 


3.3.  Optical  Hamiltonian  Regular  Equation 


.  SH  _  dp^ 
dx  dz 

dH  ^  _  dp, 
dy  dz 

dH  ^ 

^Px  dz 
dH  ^  dy^ 
dp  ^  dz 


(3.9) 


Optical  Hamiltonian  argum.ent: 

H=  -[  n-  {x,y,z)  -p\-p]V'^  (3.10) 

By  using  the  analogical  method,  the  following  useful 
conclusions  can  be  obtained:  (1)  as  is  well  known,  when  fr=h/2n 

(  h  is  Planck's  constant) - >0,  the  results  in  quantum,  m.echanics 

are  consistent  with  the  results  in  classical  m.echanics.  Upon 

analogy,  when  k’^ - >0  (short  vjavelength  approxim.ation)  ,  the 

transition  from,  wave  m.otion  optics  to  light  ray  optics  is 


32 


executed.  (2)  With  the  principle  of  indeterminacy  in  quantum, 
m.echanics , 

AxAp^^Till  (3.11) 


(ClK  and  £lp  are,  respectively,  the  coordinate  and  m.omentum, 
indeterminates )  with  analogy,  in  laser  optics  there  is  [12] 


(3.12) 

In  the  equation,  6^  and  VJgj  are,  respectively,  the  far-field 
divergence  angle  and  the  light  beam  waist  dimension  of  the  laser 
light  beam  with  mode  sequence  number  m.;  m.“  is  the  correlated 
factor  of  the  light  m.ode  (referred  to  as  the  M  -vector)  .  Only 
for  a  gaussian  light  beam,  of  ideal  fundamental  mode  {m=0)  ,  there 
is 


00^0=-^  (3.13) 

Eqs .  (3.12)  and  (3.13)  explain  that  small  light  beam  divergence 

and  small  light  beam,  waist  dimension  are  unable  to  be  realized  at 
the  same  tim.e;  this  is  the  fundamental  relationship  equation  in 
controlling  light  beam,  quality.  Eq.  (3)  m.arJ^s  the  derivation  of 
the  quantUiii-m.echanical  form.ula  system,  by  using  operators  as  the 
m.echanical  quantity.  On  analogy,  the  operator  method  can  be 
conveniently  used  to  deal  with  the  light  beam,  transmission 
transf orm.ation  with  diffraction  effect  for  development,  thus 
perform.ing  the  operator  optics. 


IV.  Operator  Optical  Method 

Operator  optics  is  expressed  by  a  method  proposed  and 
developed  by  Nazarathy  and  Sham.ir  et  al.  [15,  16].  In  the 
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expression  method,  expression  with  a  regular  operator  is  more 
complete.  The  nonregular  operator  method  is  still  under 
developm.ent  [17]  . 


4.1.  Fundam.ental  regular  operator  and  operator  form  in  Collins' 
formula 

Let  us  assum>e  that  u(x)  is  a  random  function.  From.  Eq. 
(3.6),  we  have 


xuix )  =  .Y  u(x) 


P,uix) 


-ik-' 


dujx) 

ex 


(4.1) 


and  (4.2) 

Define  the  regular  transf orm.ation  operator  T,  which  satisfies 
(for  convenience  in  writing,  the  subscript  x  is  om.itted  in  p^  in 
the  following) : 


,X:=  D.y,  -  Bp^ 
p.  =  -Cx^■^  A  p, 


(4.3) 


In  the  equation.  A,  B,  C,  and  D  from 


-Y, 


Pz 


are  the  M  elem.ents  of  the  transformation  matrix  in  the 
expression . 


It  is  easy  to  prove  that  T  satisfies 


(4.4) 


T[  M,]T[  A/,]  = 

}T[A/,]  =T[A/3]{[A/,]T[Af,]  } 


(4.5) 

(4.6) 
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(4.7) 


{T[M]}-'  =  T[A/-'] 

T[E]-1  =  °j) 

Therefore,  the  regular  operators  T  constitute  a  group.  Actually, 
T  is  part  of  the  Metaplectic  group,  isomorphous  with  the 
Symplectic  group. 

In  the  following,  the  fundamental  regular  operators  can  be 
expressed  in  terms  of  T  as 

(1)  Square  phase  operator 

2[a]=exp(-^ax^)=T[’  J]  (4.9) 

(2)  Scale  operator 

V[)?]i/(.v)=t[^  '  «(x )=)?'” (4.10) 

(3)  Fourier  transform,  operator 

r*”  n  1 

yx,exp( -/A:x,x)---  =  rj^_  J  pj  (4.11) 

V  -05 

(4)  Fennell  free-space  transm.ission  operator 

R[7]  =exp(  - =  (4.12) 

By  using  Eqs .  (4.9)  to  (4.12),  as  well  as  the  related  com.putation 

equations,  Collins'  form.ula  (Eq.  (2.10))  can  be  written  as 
(neglecting  the  nonessential  phase  factor  e  ,  hereinafter) 

(^2(-V2)=Q[  D/B]V[1/B]FQ[A/B]U,{x,)  ,  (B^O)  (4.13a) 

=Q[C/  A]W[l/ A]U,(x,)  ,iB^0)  (4., 3b) 


35 


4.2.  Applications  of  special  optical  system 


By  using  the  fundam.ental  regular  operators,  it  is  easy  to 
sim.plify  the  com.m.only-seen  optical  system.s  in  physics  and  optics, 
such  as : 

(1)  Fourier  transform,  system. 

(2)  im.aging  system. 

(3)  rem.ote  focus  system 


D  '  B  B 

^[o 

(4)  Fennell  transform  system 

T  r  ■®~|  =V[  —  ]FQ[  —  1 

L-5-'  oJ  ^  ^  5 

(5)  ideal  am.plifier 


(6)  ideal  spectrum  analysis  instrument 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


4.3.  Misalignment  optical  system. 

Let  us  assume  that  a  transf orm.ation  matrix  is  the  optical 

system,  of  (  ^  ^  .  Relative  to  the  position  and  direction  of 

CD 
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the  ideal  optical  axis,  the  misalignments  are,  respectively, 
and  Mp  (as  shown  in  Fig.  1).  the  misalignm>ent  vector  is 


(4.20) 


By  introducing  the  linear  phase  alignment  G  and  the  drift 
operator  S  [18] 

G[<5]u(x)=e'*‘^«(x)  (4.21) 

S[  m]uix)=u(x-m)  (4.22) 


the  misalignment  operator 

W[^]=W[ 


(4.23) 


The  Heisenberg-Weyl  transform  operator 

W[  /i./j  =exp  [ik(l+  Y  [  -/i,]G(  -/ij 


and  the  augmented  matrix 


A  = 


A  B  E  \ 
C  D  F  \ 
0  0  1/ 


(4.24) 


(4.25) 


Later,  by  using  the  operator  comvm.utation  equation,  the 
misalignment  transformation  operator 

H  =  W-'[  /iJT  [M]W{#i,] 


can  be  written  as  [19] 
1  ,  -^^2 


H  =  T 


0  0 
=  T[A,L] 


1 


\  1 

7  M  ,  0  ' 

\  n" 

j.y 

T 

Lvo  0  '  1  , 

).oJ 

T  1 

0  0  1 


1 

2  A*  1 i/> 

(4.27) 
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rig.  1  .A  misalignmcnl  opiic.il  >\sicni 


In  the  equation,  the  total  augmented  matrix  is 


1 


(4.28) 


the  phase  drift  factor  is 


(4.29) 


Up  to  now,  great  successes  have  been  obtained  by  using  the 
regular  operator  expression  method  to  deal  with  problems  of 
coaxial  and  misalignment  optical  systems,  physical  optics, 
Fourier  optics,  and  optical  resonators.  A  shortcoming  is  the 
fact  that  the  "aperture  factor"  A  is  used  to  express  the 
diaphragm  of  the  hard-edged  small  aperture.  This  aspect  needs 
further  study. 


4.4.  Regular  operator  expression  of  time-domain  diffraction 
integral 

Similar  to  using  the  regular  operator  method  in  the  space 
domain,  in  the  time  domain  we  also  can  define  the  regular 
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conjugate  operator  and  the  regular  transform  operator.  Then  the 
time  domain  diffraction  integration  equation  (2.12)  can  be 
written  as  the  following  equations  by  using  the  time  domain 
square  phase  operator  Q,  Fourier  transform  operator  F,  and  the 
scale  operator  V  [11] 

t/2(^2)=Q[^]V[y  ]FQ[^]t/,(t.).  (5#0)  (4.30a) 

'^2(^2)=Q[  (4.30b) 

When  using  the  time  domain  diffraction  integration  in 
studying  transmission  problems  of  chirp  pulses  in  dispersion 
medium,  it  is  not  required  to  follow  the  conventional  method  for 
multiple  transformation  of  the  time  domain  and  the  frequency 
domain.  Furthermore,  Eq.  (4.30)  written  in  terms  of  operators 
can  have  matrix  and  operator  commutation  operations  to  replace 
the  time-consiiming  integration  computations  in  the  intermediate 
steps;  therefore,  the  problem  is  greatly  simplified,  with  the 
present  form,  such  as  [18]: 

t/,(T,)=£'„exp(-2In2Tj/A;f)  (4.31) 

The  time  transformation  matrix  of  chirp  pulses  ((^Ij  is  pulse 
duration  and  Eq  is  constant)  is 


M  =  ( 


C, 


)  =  ( 


1 

0 


(4.32) 


By  using  transmission  in  a  dispersion  medium  for  explanation, 
substitute  Eq.  (4.32)  into  Eq.  (4.30a),  and  make  a  computation 
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(4.33) 


operation,  thus  obtaining 


.D,  /#,  A  A 


In  the  equation. 


A  B 
C  D 


l  +  i- 


41n2 

41n2 

(Do  A/^ 


5,  B, 


(4.34) 


Thus,  we  obtain 


exp[- 


4^(1+ ^  ) 


TPi - ]•  exp[  «>„„,( ^2)]  (4.35) 


In  the  equation. 


<Pouc(^2)  =  [  -(p"tj/(2(p"^+8Q^)]  -  y  arctg(  - 


9" 

2n 


) 


(4.36) 


Q  =  Arf/8  1n2 
<p"=P"z 


(4.37) 

(4.38) 


The  results  obtained  are  the  same  as  those  in  [20]  and  [10]. 
From  Eqs .  (4.31)  and  (4.35),  we  know  that  after  passing  through 
the  pure  dispersion  medium,  the  chirp  pulse  duration  (!^t  is 
broader,  and 


A/; 

At, 


=  (1  + 


1. _  \\n 


4Q' 


(4.39) 
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With  respect  to  the  application  examples  for  B=0  in  Eq.  (4.30b), 
these  examples  can  be  found  in  reference  [11]. 

Remark;  some  contents  of  the  paper  were  reported  in  a 
symposiiom  on  the  ICF  driver  technology  in  Beijing  in  July  1991. 
The  first  draft  of  the  paper  was  received  on  September  12,  1991; 
the  final  revised  draft  was  received  for  publication  on  June  15, 
1992. 
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